ADb58S45 


a. 


:  ^ 


■hi 

/ 

EFFECTS  CF  ABLATION 

fly 

ON  THE  \ 

PITCHING  MOMENT  DERI  VAT 

IVES 

OF 

|  CONES  IN  HYPERSONIC  FLOW 

V-- 


.  V  — .V 


'  ■.  r. 


):  \ 


r , 


\ 


_,  _  V 

NP*  . 


./  . 


*  ^ 


•/ 


ups? 

AVAILABLE  COPY 


FlUIDYN!  ENGINEiAING  CORPORATION 


■^5  > 


PITCHING  MOMENT  DERIVATIVES 
OF 

CONES  IN  HYPERSONIC  FLOW 


by 

Shukry  K.  Ibrahim 

FltiiDyne  Engineering  Corporation 
5900  Olson  Memorial  Highway 
Minneapolis,  Minnesota  55422 
Report  0414-4-4-67 


This  work  is  being  supported  by  the 

Advanced  Research  Projects  Agency 
Ballistic  Missile  Defense  Office 
Pentration  Aids  Branch 


and  technically  administered  by  the 

Fluid  Dynamics  Branch 
Office  of  Naval  Research 


Order  Number  ARPA  No.  576 
Contract  Nonr-4624 ( 00 } 


D  D  C 

r~3  r — -  --  •  -  '“irp 

Dv 

SEP  1  «  1967 

Ulta~ _ 1Lj  U  i-lj 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government. 


Distribution  of  this  document  is  unlimited 


C\ 

N 


FluiDyne  Engineering  Corporation 


FOREWORD 

The  work  described  in  this  report  was  carried  out 
as  part  of  a  program  to  investigate  the  effects  of  mass 
transfer  on  the  stability  of  re-entry  type  vehicles.  In 
addition  to  the  theoretical  analysis  reported  here,  an 
experimental  investigation  was  conducted  in  the  FluiDyne 
Hypersonic  Wind  Tunnel  at  a  nominal  Mach  number  of  11. 

It  involved  free  oscillation  tests  with  ablating  and 
blowing  models  using  a  specially-designed,  sting-mounted 
air  bearing  with  the  necessary  control  and  recording 
apparatus.  The  overall  program  was  supported  by  the 
Advanced  Research  Projects  Agency,  Ballistic  Missile 
Defense  Office,  Penetration  Aids  Branch,  and  was  tech¬ 
nically  administered  by  the  Fluid  Dynamics  Branch  of 
the  Office  of  Naval  Research.  The  author  wishes  to 
acknowledge  the  essential  contribution  of  Dr.  James  S. 
Holdhusen,  the  co-author  of  an  earlier  report  on  the 
subject  which  laid  the  foundation  for  the  present  work. 
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ABSTRACT 

This  work  presents  a  method  of  calculating  the  static 

<■ — 

and  dynamic  pitching  moment  derivatives  at  hypersonic  speeds 
for  pointed  or  slightly  blunted  cones  performing  small- 
amplitude,  s i ngle-aegree-of-freedom  oscillations  about  an 
arbitrary  pivotal  axis,  taking  into  account  the  effects  of 
ablation  or  blowing  over  all  or  specified  portions  of  the 
cone  surface.  The  study  is  motivated  by  the  need  for  de¬ 
tailed  understanding  of  the  large  destabilizing  changes 
in  the  pitching  oscillations  sometimes  observed  on  re-entry 
vehicle  test  flights.by  developing  analytic  expressions 
which  describe  the  role  and  relative  importance  of  the 
relevant  parameters. 

An  earl>  report  on  the  subject^  with  a  limited  dis¬ 
tribution  was  presented  in  April  1966.  Significant  im¬ 
provements  have  since  been  made  to  the  earlier  paper  in 
that  the  effects  of  skin  friction,  as  modified  by  ablation 
or  blowing,  on  the  static  and  dynamic  pitching  moment  de¬ 
rivatives  are  now  included.  Other  refinements  include  the 
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choice  of  a  body  coordinate  system  centered  at  the  cone 
apex,  or  virtual  apex  for  blunted  cones,  leading  to  a 
more  unified  approach  to  the  effects  of  ncse  bluntness 
and  the  relaxation  of  certain  geometric  restrictions. 

The  analysis  is  restricted  to  considerations  of  rigid 
body  motions  at  low  values  of  the  reduced  frequency  param¬ 
eter,  ,  and  quasi-steady  conditions  are  postulated. 

00 

The  formulation  adopted  is  sufficiently  general  and  per¬ 
mits  a  great  deal  of  latitude  in  the  choice  of  parameters 
that  apply  to  specific  cases.  The  primary  independent 
quantities,  other  than  the  cone  angle,  nose  bluntness  and 
center  of  gravity  location  are  the  distribution.  Intensity 
and  phase-shift  of  the  ablation  or  blowing  rate.  In  the 
analysis,  the  ablation  gases  leaving  the  body  are  assumed 
to  produce  the  following  two-fold  effect: 

1)  a  deflection  of  the  external  flow  which  produces 
a  chruge  in  the  local  static  pressure  coefficient 
readily  calculated  by  the  Newtonian  flow  theory  ar.d 

2)  a  change  in  the  local  skin-friction  coefficient. 

These  inviscid  and  viscous  effects  are  treated  separately, 
and  their  individual  contributions  to  the  static  and  dy¬ 
namic  derivatives  are  calculated  so  that  they  may  be  com¬ 
pared  with,  and  later  added  to,  the  static  and  dynamic 
derivatives  without  ablation.  The  resulting  expressions 
conform  to  tie  observed  sign  and  order  of  magnitude  of 
experimental  observations.  Under  certain  restrictions 
clearly  delineated,  a  very  simple  relation,  involving 
only  the  reduced  frequency  parameter  and  the  phase-shift 
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angle,  is  shown  to  exist  between  the  dynamic  and  static 
moment  derivatives  due  to  ablation. 

Illustrative  examples  showing  detailed  calculation 
procedures  and  typical  results  for  three  diffe'ent  abia 
tion  configurations  under  both  laminar  and  turbulent 
boundary  layer  conditions  are  presented. 


■) 

f 
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I  a§ 


1.  INTRODUCTION  AND  LITERATURE  REVIEW 

It  is  known  that  ablation  may  produce  order-of- 
magnitude  changes  ir.  the  dynamic  pitching  derivatives 
of  cones  in  hypersonic  flow.  Such  changes  have  been 
observed  both  in  free  flight  and  in  wind  tunnel  testsz 
with  ablation.  In  examining  the  dynamic  stability  of  an 
ablating  hypersonic  vehicle,  one  is  faced  with  very  com¬ 
plex  aero-thermo-chemical  phenomena  involving  unsteady 
heat  and  mass  transfer  processes,  multi-component  boundary 
layer  effects,  different  mechanisms  of  decomposition  of 
the  ablating  material,  and  finally,  changes  in  shape  and 
mass  distribution  resulting  from  the  removal  of  ablated 
mater i a  1 . 

Any  analysis  <•?  reasonab  1y  wi  de  application  entails, 
of  necessity,  a  number  of  simplifications  dictated  e'ther 
by  the  underlying  assumptions  of  the  theories  involved  or 
by  the  desire  to  reduce  the  necessary  computational  work. 
The  latter  restrictions  may  be  relaxed  at  the  expense  of 
additional  computational  labor  but  the  former  are  intrin¬ 
sic  to  the  nature  of  the  theoretical  approach  adopted. 

From  the  fundamental  viewpoint,  the  present  analysis 
uses  the  classical  concepts  of  dynamic  stability  theory 
such  as  those  of  stability  derivatives  and  the  principle 
of  superposition  implying  the  assumption  of  linearity. 

It  also  rests  on  the  assumptions  of  the  Newtonian  flow 
theory  and  of  quasi-steady  flow  conditions. 

From  the  practical  computational  viewpoint,  a  num¬ 
ber  of  simplifying  assumptions  are  made  in  the  interest 


\ 
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of  clarity  and  simplicity.  Only  the  primary  effects  tak¬ 
ing  place  in  the  ablating  region  are  considered  and  second¬ 
ary  effec+s  as  might  ue  induced  on  the  body  upstream  or 
downstream  of  the  ablating  region  are  ignored.  The  con¬ 
vective  heat  transfer  rate  is  assumed  to  be  dominant  so 
that  radiation  effects  may  be  ignored.  Possible  changes 
in  shape  or  mass  distribution  resulting  from  ablation  are 
also  ignored,  and  the  dynamic  stability  effects  of  ablation 
on  the  small,  spherical  nose  cap  are  neglected.  !n  prin¬ 
ciple,  there  are  no  difficulties  in  relaxing  anv  or  all  of 
these  simplifying  assumptions. 

The  open  literature  treating  spscifically  the  complex 
interactions  between  the  phenomena  of  ablation  and  the  ve¬ 
hicle  dynamic  stability  characteristics  is  quite  limited 
and  of  very  recent  vintage  .  By  contrast,  the  literature 
on  some  of  the  underlying  problem  areas  is  quite  axtensive 

and  only  a  representative  number  of  references  in  the  areas 

7-20 

of  high  speed  heat  and  mass  transfer  ,  stability  charac- 

2  1-25 

teristics  at  hypersonic  speeds  ,  ablation  processes 

2  6-  3  7 

and  material  properties  are  l:sted  here. 
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2.  PITCHING  MOMENT  DERIVATIVES  FOR  NONABLATING  CONES 


2.1  Geometric  Considerations 


The  cone  geometry  is  illustrated  in  Figure  1.  The 
origin  of  the  cylindrical  coordinate  system  (x,  r,  cp) 
is  taken  at  the  apex  of  the  cone.  For  blunted  cones, 
the  origin  is  taken  at  the  virtual  apex  as  shown. 

The  significance  of  the  symbols  L,  LQ  ,  ,  xQ  , 

r*N  >  rg  »  and  6  is  illustrated  in  Figure  1.  The  nose 
bluntness  ratio  is  defined  by: 


The  following  relatiins  apply  to  a  right  circular  cone 
with  a  spherical  nose  cap: 


r 

L 


=  xtar,6  ,  rg  =  Ltan0 

rN(l-sin^0) 

1  =  STnS  '  5  = 


,  d  =  2l,tan0 


(1) 


The  local  pressure  coefficient  for  hypersonic  flow 

for  a  cone  is  based  on  Newton's  simple  impact  theory 

{valid  for  M  -*  »)  to  be  modified  by  a  corraction  factor 

38 

^  given  in  Truitt  by  the  expression 
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r\  = 


l+]  gg) 

4(y-D  M  2A2' 


(2) 


where 

and 


Y  is  the  ratio  of  specific  heats 

hi  is  the  free  stream  Mach  number 
00 

A  represents  the  angle  between  the  tangent 
to  the  local  surface  and  the  free  stream 
d i rect i on. 


Under  these  conditions,  the  local  pressure  coefficient  is 
liven  by: 


C  (o)  =  2A2t| 


(3) 


The  local  deflection  angle,  A  ,  Is  expressed  in  terms 

of  the  cone  angle  0  and  the  local  effective  angle  of  , 

attack,  a  ,  composed  of  the  instantaneous  angle  of  at- 
6 

tack  and  the  angle  induced  by  the  P'tching  velocity  q  : 


a 

e 


a  + 


q(x-xQ) 


Introducing  the  dimensionless  quantities: 

"x=rj"xo=r'[i=r 

then  %  =  a  +  =  a  +  Kr  (tSThD  *  {4) 

GO  QO 

Following  the  earlier  analysis^.,  we  use 

A  =  0  -  a  sin  (p  (5) 

e 
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and  obtain  for  the  pressure  coefficient  with  no  ablation 
using  Eqs.  (3)  and  (5): 

Cp  ( o )  =  2r\(  0  -  sin  cp ) 2  .  (6) 

2.3  Pitching  Moment  Coefficient  and  Derivatives  for 
Nonab  1  at i ng  Cone  Surface 

The  differential  pitching  moment  coefficient  about 
the  pivotal  axis  located  at  xQ  for  a  nonablating  conical 
element  at  x  is  given  by 


7  rcp=n/2  _ 

dCM(o'cl  =^If=.tI/2[r  Sin0  + 


(x-xo)cos8jCp(o)sincp  dS 


(7) 


where  the  cone  base  area  AD  and  base  diameter  d  are 

B  2 

used  for  reference.  Noting  that  Ag  =  nrg  ,  d  =  2rg 

and  that  the  element  of  area  dS  on  the  cone  surface  is 

gi  ven  by  dS  =  >  then  by  substitution 

cos  8 

from  Eq.  (1)  and  introducing  the  dimensionless  quantities 
x  ,  xQ  ,  and  ,  there  is  obtained  from  Eq.  (7)* 


dC^(o,c  i 


- T“ 

TTt an  8 


»*/ 2  _  t  9  _  _ 

i  (x^sec  8  -  xx  )  C  { o ) 
-tt/2  0  p 


s  i  nip  dip  dx 


The  moment  coefficient  C^(o,c)  ,  contributed  by  the 
nonab'ating  cone  surface,  requires  an  integration  with 
respect  to  x  from  the  shoulder  of  the  cone  near  the 
nose,  i.e.,  x  =  to  the  :one  base,  x  =  1 
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,1  oTt/2 


C.,(o,c)  - - f  [  (x2sec20  -  xx  )  C„(o)  s  i  ncp  dtp  dx  (8) 

M  ntan20  Jt,  J -tt/2  0  p 


Thus,  the  moment  coefficient  depends  on  the  geometrical 
quantities  0  ,  L|  ,  and  xQ  and  on  the  distribution  of 
the  pressure  coefficient  with  respect  to  the  variables 
cp  and  x  .  The  nose  bluntness  ratio,  §  ,  does  not  ap¬ 
pear  explicitly  in  Eq.  (8)  but  its  effect  is  incorporated 
in  the  lower  limit  of  integration,  . 

It  should  be  pointed  out  that,  except  for  pointed 
cones,  the  numerical  values  of  xq  ,  as  defined  in  the 
present  and  in  the  previous  report^  will  be  different 
for  the  same  geometrical  configuration  and  axial  loca¬ 
tion  so  that  the  calculated  moment  coefficients  and  their 
derivatives  cannot  be  strictly  compared  with  xq  as  the 
independent  variable. 


By  substitution  from  Eq,  (6)  into  Eq.  (8),  there 
foil ows : 


C.,(o,c)= — 2tU f  f  {x2sec20  -  xx  )  (  0-a  s  i  ncp ) ^ s  i  ncp  dtp  dx  (9) 

M  ntan20  JL1‘-tt/2  0  e 


The  expansion  of  Eq.  (9)  yield  terms  involving  s  i  nqp  and 

3 

its  powers  up  to  sin  cp  •  Since  the  integrals  of  the  odd 
powers  of  si  ncp  with  limits  cp  ~  -tt/2  and  +n/2  will 
vanish  and  since 

r"72  . 2  n 

sin  cp  dep  =  "k 

-n/2  2 
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Then,  after  integration  of  Eq.  (9)  with  respect  to  cp 

and  substitution  for  cc  from  Eq.  (4),  there  is  obtained 

6 

for  the  moment  coefficient  without  ablatio.) 

cMlo,c)  =  -  «2lL|-  aj]  (x2sec20  -  xxQ)dx  - 
tan  0  L  ^ 

2tL|-  {x2sec20  -  xxQ )  Cx-xQ )dx  .  (10) 

t  an  0  ®  L  ^ 


After  carrying  out  the  integrations  with  respect  to  x  , 
taking  appropriate  derivatives  and  using  the  abbreviations: 


-  1. 


C,  -  M 


sec2  0 


(11) 


and 


l-I 


C2  =  (  )sec2e 


O 

sec^e 


+ 


2 

0 


(12) 


"here  is  obtained  for  the  pitching  moment  derivatives: 


CM  (o,c) 
a 

and 


-tof-c, 

tan^  0  1 


CM  (o,c) 


2  ri9 
tanJ0 


( 13) 


(14) 


It  is  observed  that,  for  a  given  problem,  and  C2  are 

Constant  coefficients,  specified  by  purely  geometric  fac¬ 
tors,  namely,  the  cone  angle  0  ,  the  pivotal  axis  location 
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xQ  and  the  nose  bluntness  which  controls  Lj  .  Cj  is 
linear  and  C2  is  quadratic  in  xq  .  For  the  special 
case  of  a  pointed  cone,  Ej  =  0  and  the  coefficients 
and  C2  become,  respectively: 

Cl,=  3  sec  9  "  2  *0 

and 

C9'=  4  sec20  -  4  (soc20  +  1  )x  +  4  x  2 

L  4  3  0  L  0 

Figures  2  and  3  illustrate  the  variation  of  C|  and  C2 
with  pivotal  axis  location,  xQ  ,  for  a  10°  half-angle 
cone  for  different  values  of  the  bluntness  ratio,  §  . 
Figures  4  and  5  present  the  nonablating  static  and  dynamic 
pitching  moment  derivatives  produced  by  the  conical  sur¬ 
face  and  given  by  Eqs.  (13)  and  (14). 

2.4  Pitching  Moment  Derivatives  of  Spherical  Nose  Cap 

The  pitching  moment  csrivatives  contributed  by  the 
spherical  nose  cap  are  readily  calculated  by  means  of  the 
Newtonian  Impact  Theory.  From  Eq.  (9)  of  Reference  39 
but  with  the  cone  base  a^ea  and  base  diameter  as  refer¬ 
ence  area  and  reference  length,  respectively,  the  pitch¬ 
ing  moment  derivatives  about  xQ  ,  contributed  by  the 
spherical  nose  are  expressed  by 

CMa(o'N)  =  V?U  {*o  ‘  lo  '  5  tan6) 


(15) 
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CM  (o,N)  = 


§2cos40 

- 2 — 

tan  0 


§  t  a  n  0 ) 2 


(16) 


These  contributions  to  the  static  and  dynamic  pitching 
moment  derivatives  for  a  bluntness  ratio  §  =  .2  are 
illustrated  in  Figures  4  and  5. 


2 . 5  Pitching  Moment  Derivatives  of  Blunted  Cones  Without 
Ab lat i on 


The  resultant  static  and  dynamic  pitching  moment 
derivatives  for  the  nonablating  spherically  blunted  cone 
are  calculated  from  Eqs.  (13),  (14),  (15),  and  (16)  by 
addition,  so  that 


Cu  (o)  =  Cu  (o,c)  +  Cu  (o,N )=  -  -232-  C,  + 
Ma  Ma  Ma  tan  8  1 

<*o  -  to  -  5 


and 


Cu  (o)  =  CM  ( o , c )  +  CtJ  (o,N )=  -  -22$-  C,  - 


M. 


3  ^2 

tan'3©  z 


lLqos4i  (x  -  t  -  §  tan 0 )2 
tan2  0  0  0 


(17) 


(  18) 


The  calculated  resultant  values  for  a  10°  half-angle  cone 
with  a  bluntness  ratio  §  =  .2  are  illustrated  in  Figures 
4  and  5. 
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3.  CHANGES  IN  PITCHING  MOMENT  DERIVATIVES  DUE  TO  ABLATION 
3.1  General  Considerations 


The  general  approach  will  be  to  calculate  the  convec¬ 
tive  heat  transfer  rate  on  the  local  element  of  conical 
surface  by  suitable  transformation  of  flat  plate  data  at 
angle  of  attack,  choosing  the  latter  to  correspond  with 
the  local  flow  inclination.  The  loca1  rate  of  mass  trans¬ 
fer  of  ablation  gases  into  the  boundary  layer  may  then  be 
calculated.  Either  one  of  the  following  procedures  may 
be  adopted. 


or 


A.  To  calculate  the  heat  transfer  rates  by  means 
of  accepted  calculation  procedures  assuming  no  abla¬ 
tion  and  then  apply  the  concept  of  an  "effective 
heat  of  ablation,"  H^^  ,  to  account  for  the  thermal 
cooling  properties  of  the  ablator  and  ambient  flow 
cond i t i ons. 


B.  To  calculate  the  heat  transfer  rates  taking 
into  account  the  alleviation  resulting  from  the 
ablating  process  and  then  to  determine  the  local 
rate  of  mass  transfer  from  the  known  thermal  pro¬ 
perties  of  the  ablating  surface  mater  ial. 


The  first  method  was  adopted  and  led  to  a  simple  analysis. 
It  should  be  stressed  that  the  effective  heat  of  ablation 
is  not  a  material  property  of  the  ablator  exclusively  but 
depends  on  the  flow  characteristics  as  well. 
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In  turn,  the  local  mass  release  at  the  surface  is 
assumed  to  produce  a  two-fold  effect: 

1)  a  local  incremental  deflection  in  the  flow 
which  contributes  an  increment  to  the  local  pressure 
coefficient,  designated  by  Cp(0)  ,  and  readily 
calculated  by  a  simple  Newtonian  analysis, 

and 

2)  a  change  in  the  local  boundary  layer  velocity 
profile,  more  specifically  a  reduction  of  the  nor¬ 
mal  derivative  of  the  velocity  at  the  wall  surface 
and  a  consequent  reduction  of  the  local  skin  fric¬ 
tion  coefficient. 

Vhe  inviscid  and  viscous  effects  are  treated  separately 
and  their  respective  contributions  to  the  pitching  moment 
coefficient,  designated  by  C^O)  and  CM  ( f )  are  cal¬ 
culated  independently  and  later  added  to  the  pitching 
moment  coefficient  without  ablation,  ( o )  ,  to  yield 
the  resultant  pitching  moment  . 

3.2  He  :t  Transfer  Rates  on  Cone  Surface  Element 


Simplified  expressions  for  the  rate  of  convective 
heat  transfer  to  a  flat  plate  at  angle  of  attack,  a  , 
under  re-entry  coi  Jons,  have  been  presented  in  Refer¬ 
ence  40.  For  the  laminar  boundary  layer  case,  the  follow- 
.ng  expression  is  given: 


.  3068  _ sin  a  ( cos  a)  ‘  ^ 

/  Y_  i  \  •  2  5  "Z  h  2  •  .25 

J  .  1 8  6  +  .5  -pj—  +  .  3 1 4  ( s  i  not ) 


(  19) 
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where 

=  convective  heating  rate  for  1~  inar  flow 
(BTU/f t2sec ) 

x  =  coordinate  distance  along  the  surface  (ft) 

a  =  altitude  density  ratio  ,  £» 

PS.L. 

=  free  stream  velocity  (ft/sec) 
h  =  enthalpy 

Y  =  ratio  of  specific  heats 

w,s  =  subscripts  to  designate  the  wall  and  stag¬ 
nation  values,  respectively. 

For  small  values  of  a  ,  such  that 

sin  a  »  a  and  cos  a  &  1 


it  is  shown  in  Appendix  B  that  one  may  use  the  approxi- 
mat i on 


.5 

s i n  a  (cos  a) 

r  FT  ‘  5^75 

•  1 86  +  .5  -jr-  +  . 3 14 ( s  i  n  ot)  J 


(20) 


The  heat  transfer  rate,  ^ ( F , P . )  ,  on  the  flat  plate  may 
be  converted  to  that  on  the  cone,  i)(c)  at  the  same  dis¬ 
tance  x  by  means  of  the  Mangier  transformation.  For 
the  laminar  boundary  layer  case: 

qt(c)  =  J5  ^(F.P.  )  .  (21  ) 

It  should  be  pointed  out  that,  for  the  cone,  the  flow 
conditions  must  be  those  at  the  edge  of  the  boundary 
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layer,  designated  by  subscript  e  ,  and  not  the  free 
stream  conditions  indicated  by  subscript  »  .  The  dis¬ 
tance,  s  ,  along  the  conical  surface  for  a  point  with 
coordinate  x  is  given  by  s  =  x/cos0  =  Lx/cos0  .  The 
instantaneous  local  inclination,  A  ,  of  the  flow  which 
corresponds  to  the  angle  a  of  the  flat  place  is  given 
by  Eq.  (5). 

By  virtue  of  Eqs.  (19),  (20),  (21),  and  (6),  there 
is  obtained  for  the  local  heat  transfer  rate  on  the  coni 
cal  surface  under  laminar  flow  conditions: 


h.\ 

V 


.3068 

"23 


L_ ' 5(cos0) *  5 

7 - h“TT73 

[•186+  -5^] 


Ax 


-.5 

c^d-x) 


(22) 


where  the  coefficient  C.  is  given  by 


(cos  0) 


.5 


(¥r5  [. 


-_2_ 

186  +  .5 


(23) 


It  is  determined  by  the  geometric  configuration  and  the 
flow  character i st i cs . 


The  local  heat  transfer  rate  on  the  conical  surface 
under  turbulent  boundary  layer  conditions  may  be  calculated 
along  similar  lines.  Reference  40  gives  the  simplified 
heating  equation  for  the  flat  plate  at  angle  of  attack  for 
the  turbulent  case  as 
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1.6  .8 

( s i n  a)  ( cos  a) 

7 - k - 

[.197+  .5  -jp  +  . 3 03 ( s i  na)  1 


(24) 


It  is  shown  in  Appendix  B  that,  for  small  values  of  a  , 
one  may  use  the  approximation: 


1.6  .8 

(sin  a)  (cos  a) 

F 


1.6 


— TT7 

[.  197  +  .  3  03  ( s  i  n  a )  J 


.  197  + 


K  t  ~7 

5ttJ 

s 


(25) 


1 

■I 

* 

3 


The  Mangier  relationship  for  this  case  yields 
qt(c)  =  1.14  $ t  (F.P. ) 


(26) 


Following  the  same  procedure  as  before,  there  is  obtained 
for  the  local  heat  transfer  rate  on  the  conical  surface 
under  turbulent  boundary  layer  conditions: 


*t<‘>  -  -8(fc5of('  -  r-',2|w6lhb  *K6  *'-2 

S(-V)  [.197  +  .5-^] 

o 


Ct  A 1 ’ 6 ( Lx ) 


-.2 


(27) 


where 


Ct  =  1.  14  o 


.6 


/  U  N3/ 

(lo§o)  ( 


h  N 

-  -2*0 
h  _/ 


3.711  (cos  Q)'2 


(28) 
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In  o^der  to  unify  the  analysis  of  the  laminar  and 
turbulent  cases  and  generalise  the  approach  (to  encompass 
other  ablation  or  blowing  configurations),  it  is  conveni¬ 
ent  to  adopt  the  general  formulation 

q(c)  =  C  Am ( Lx ) P  (29) 

where  the  coefficient  C  represents  an  intensity  factor 
the  exponent  m  indicates  the  angular  dependence 
and  the  exponent  p  denotes  the  x-wise  variation  of 
the  heat  transfer  rate. 

In  general  m  is  positive  but  p  is  negative;  the  special 
case  p  =  0  implies  no  direct  x-wi se  dependence  of  the 
heat  transfer  along  che  surface. 

For  the  laminar  case: 

C  =  ,  m  =  1  .  and  p  =  -.5 

and  for  the  turbulent  case 

C  =  Cj.  ,  m  =  1.6  ,  and  p  =  -.2 

3  •  3  Mass  Transfer  and  Flow  Deflectiot  •  _ju  c  to  Ab  lat  i  on 


We  introduce  the  dimensionless  mass  injection  param¬ 
eters: 


and 


F  = 


F_  = 


P  v 
rw  w 

p  b 
re  e 


P  v 
rw  w 

pTT 

*  00  ea 


(30) 


(31  ) 
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where  v  represents  the  normal  velocity 

U  represents  the  tangential  velocity 
p  represents  the  density 


and  the  subscripts  w  ,  e  ,  and  08  indicate  conditions  a‘ 
the  wall,  edge  of  the  boundary  layer,  and  the  free  stream, 
respectively.  Applying  the  concept  of  the  effective  heat 
of  ablation,  H  ,  we  write 


qk)  =  x  Heff  =  (Pwvw)  x  Heff 


By  virtue  of  Eq.  (29),  there  is  obtained 


=  C  Am(Lx) 


m 


w 

Therefore 

and 


H 


eff 


-  C  ( Lx ) 

Fe  "  Heff  W 


»m  f  ■  -  * p 


C  A  ( Lx  ) 

»  ~  K'  (p  Um[ 

eff  00  09 


It  is  to  be  noted  that 


(32) 


(33) 


P  U 

*  GO  GO 

Fe  =  pir  F-  *  bF- 
re  e 


(34) 


pJJ, 

where  b  is  used  t>"  designate  the  quantify  -y  which 

pe  e 

may  be  considered  a  constant  quantity  for  any  given  prob 
lem  being  a  function  of  the  free  stream  Mach  number,  M^ 
and  the  cone  angle. 
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The  additional  local  flow  deflection  due  to  abla¬ 
tion  mass  transfer  is  designated  by  P  .  From  the  wot  x 
of  Mann*^,  Thyson  and  Schurmann^  and  others,  one  may 
write  to  a  satisfactory  approximation 


F  + 
e  dx 


(35) 


where  6*  designates  the  displacement  thickness  of  the 

boundary  layer  with  ablation  or  blowing.  It  is  fair  to 

1 9 

assume,  following  Oanberg,  Winckler  and  Chang  and  also 

6 

Ericsson  and  Reding  that  the  rate  of  increase  cf  the 
displacement  thickness  with  distance  is  linearly  related 
to  the  mass  injection  parameter  so  that 


P  =  a  Fe  .  (36) 

Experimental  evidence  regarding  the  proportionality  fac¬ 
tor,  a  ,  indicates  its  dependence  on  the  flow  characteris¬ 
tics,  more  specifically  on  the  wall  temperature  and  flow 

Mach  number  and  may  range  between  6  and  16.^  For  example, 
18 

Oanberg  obtained  a  value  of  8.5  for  turbulent  flow  on  a 

43 

flat  plate  at  M  =  6.7  in  air  and  Studerus  obtained  a 
value  of  6.2  for  turbulent  flow  over  a  cone  in  air  at 
Mach  10. 


By  virtue  of  Eqs.  (33),  (34),  and  (36),  there  is  ob- 
ta  i  ned : 


0  = 


a  b  C  Am { Lx ) 


eff 


Tpjrr 


(37) 


^  The  author  is  indebted  to  Messrs.  M.  L.  Roberts  and 
J.  B.  Arnaiz  for  information  and  references  about  the 
proportionality  factor,  a  ,  and  its  rar.ge  of  variation. 
This  information  was  given  in  prepared  comments  on  the 
author's  earlier  paper. 
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The  quantity,  K,  = -n - -it-  ,  is  a  constant  for  a 

1  neffp®» 

particular  configuration  and  with  this  abbreviation,  the 
expression  for  the  deflection  angle  due  to  ablation  be¬ 
comes 


P  =  Am  xP  .  (38) 

One  final  consideration  remains  to  be  introduced,  namely, 
that  of  possible  phase-shifts  between  the  instantaneous 
local  flew  character i s t  i  cs  at  the  surface  and  the  conse¬ 
quent  mass  release  due  to  ablation;  there  may  also  be 
phase-shifts  between  the  local  mass  release  and  the  ensu¬ 
ing  flow  deflection.  In  this  analysis,  the  combined  ef¬ 
fects  of  thermal  diffusivity  and  viscous  lag  are  assumed 
to  produce  a  time  shift.  At  ,  considered  to  be  constant 
ovo-  the  full  oscillation  cycle  so  that  the  angular  phase- 
shift,  designated  by  A  will  De  given  by 

A  -  — radians  (39) 

where  T  is  the  period  of  oscillation. 

The  determination  of  the  actual  phase-shift  for  a  partic¬ 
ular  problem  requires  a  separate  analysis  which  may  involve 
consideration  of  unsteady  boundary  layer  and  temperature 
change,  thermal  diffusivity  properties,  etc.  Such  an 
analysis  goes  beyond  the  scope  of  this  work.  Here,  aver¬ 
age  phase-shifts  involving  either  a  lag  or  a  lead  with 
respect  to  the  angle  of  oscillation  are  assumed  and  their 
possible  effects  on  the  stability  derivatives  are  examined. 

The  expression  for  the  angle  8  taking  account  of 
a  phase-shift  angle  with  respect  to  the  oscillation 
angle  a  becomes: 
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0  =  K1  Am  xP  e .  (40) 

For  an  actual  re-entry  vehicle  following  a  specified 
trajectory,  the  vehicle  speed  as  well  as  the  ambient 
atmospheric  conditions,  e.g.,  pressure,  density,  etc., 
will  continuously  change.  However,  it  is  convenient 
and  quite  adequate  for  the  purpose  of  calculating  the 
dynamic  stability  characteristics  to  assume  that  the 
re-entry  cone  is  moving  at  some  representative  velocity 
and  altitude  conditions  while  undergoing  sma  1 1-ampl  i  t  jde, 
low-frequency,  essentially  sinusoidal  oscillations  about 
the  zero  angle.  Such  are  also  the  conditions  simulated 
in  the  free  oscillation  wind  tunnel  tests.  In  terms  of 
an  initial  displacement  angle  clq  and  oscillation  fre¬ 
quency  io  ,  the  angle  a  and  pitching  velocity  q  a'-e 
expressed  by: 

a  =  a  e'“)t  ,  q  =  fit  =  iiua  eiuJt  =  iioa  .  (41) 

o  o 

Eq.  (6)  for  the  pressure  coefficient  without  ablation  may 
be  written  in  the  abbreviated  form: 

C  (o)  =  2q 02y2  (42) 

o  sincp 

where  y  =  1 - g - 

Noting  that  A  =  0y  ,  Eq.  (40)  becomes 

0  =  K1  0m  ym  xP  e'A  .  (43) 

The  Newtonian  pressure  coefficient  with  ablation  nay  be 
expressed  by 

Cp  =  2r,  (0y  +  0)2  .  (44) 
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Therefore,  the  pressure  coefficient  contributed  by  ab  1  a- 
tion,  CpO)  ,  may  be  written  as 


VP1  =  CP  -  Vo) 


By  virtue  of  Eqs.  (42),  (43),  and  (44),  one  obtains 


,m+l  m+l-p  i'A 


2 ft2m. 2m-2p  2iA 


CO)  =  4ti  K |  9 '  ym  x  e‘A  +  dr\  K/sT’x  e 


(45) 


aesincp 

For  I  — g -  |  <  1  ,  the  binomial  series  expansion  yields: 


m 


a 


s  i  ncp 


m 


)  -L  <-'>"(*)  (■ 

n=0 


sincp.n 


e—1) 

0  ) 


When  this  is  applied  to  Eq.  (45)  there  is  obtained  for 
the  pressure  coefficient  contributed  by  ablation 


CpO)  =  4r|  K,  0m+1  xPeiX  Y  (-1)"  (m+])  (-e-g--CP)n 

n=0 


+  2n  K,2e2m  -x2°  e2i*  l  1-1)"  (2m)  (V^)"  . 


(46) 


n=0 


3.4  Pitching  Moment  Due  to  Ablation  Mass  Transfer 


The  expression  for  the  pitching  moment  coefficient 
contributed  by  ablation  takes  the  same  form  as  Eq.  (8) 
but  with  CO)  in  the  place  of  C  ( o )  and  with  the 

r  —  r 

appropriate  limits  for  the  x- i ntegrat i on.  In  other 
words 
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,X2  ptt/2  _2  o 

CU(P)  - - w—  (x  sec  9-x  x  )C  (0)  sincp  dcp  dx  .  (47) 

M  ntari  9  Jx,  J-tt/2  0  P 


Consider  first  the  integration  with  respect  to  cp  : 


jt/2  n+ 1 

(s  i  ntp)  dcp  = 

-tt/2 


V-A  ,P.°s.T  |n/2  +  n  r"/2  d 

^+1  '.n/2  n+1  J.„/2  9 


The  first  term  on  the  right-hand  side  vanishes  for  the 
integration  limits  -tt/2  and  +tt/2  and  the  second  term 
vanishes  only  for  even  values  of  n  ;  for  odd  values, 
one  gets: 

TT  3  5 

j  for  n  =  1  ;  -^tt  for  n  =  3  and  -y^n  for  n  =  5 


After  summation  to  n  =  5  and  integration  with  respect 
to  cp  ,  one  obtains: 

ptt/2  m+l_p  iAr  a 

1^  =  J  ^  Cp  (P)  s  i  ncp  dcp  =  4r|  K  |  9  x  e  (m+1)  y  -g - 

/  m+ 1  \  3  >  e  ^  im+li  5  >  e  i  ^  ~]  . 

(  3  1  S”  (T'  "  {  5  1  TS*  (T_)  '"}  + 


K2e2m  x 


;2P.2,X[-2.S  £) 


0  Q  a  3 

<3 


0  t-  a  5 

l5m)  TS'V1  ••  ‘J 
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3.4. 1  Laminar  Boundary  Layer 

For  the  laminar  case,  m  =  1  and  only  the  binomial 
coefficients  (m  +  1)  and  2m  will  remain  in  Eq.  (48)  so 
that 


p  i  X 


2  _2p  2iX 


I  R  =  -4t|tt  Krx  e  ae  -  2nrr  KR  9  x  e  ag 


where  the  subscript  l  has  been  used  to  indicate  laminar 
values.  By  substitution  from  Eqs,  (4)  and  (49)  into  Eq. 

(47)  there  is  obtained: 

4ti  K , .  0e 1  ^  rX2/_2+p  2  _  l+p_  w  n. 

CM(P,^)  =-  — - — 2" -  I-  (x  sec  0_x  xo)(a+  2“U?a ( x- x o ) ) d x 

tan  8  x. 


2*1  ^ 0e  f  2/_2+2p  2  _  l+2p_  w  _  \  _ 

- -IV -  J_  (x  sec  0-x  xoXa+  7ufaiTC(x~xo,)dx  *{5° 

tan  8  x, 


3+p  _  3+p 

/  X<i  —  x  i  •> 


2+p  _  2+p 


"  X 


L.‘  C3  ■  (  2  3  +  p 1  W*  - 


C  4=(- 


_  3+2 p  _  3+2 p 
'x2  "  X1  \  2 


2+2 p  _  2+2 p 


_^_)sec2e  - 


4+p  4+p  3+p  3+p 

^5  -  (  W-29  '  (W)'-2^»x0  ♦ 


2+p  2+p 


fX2  "  X1  \z 

\ — tt~p — r< 
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and  Cg  = 


'o  4+2p 
-  x  , 


3+2 p  _  3+2 p 

,x„  -  x 


2p 


)sec20  ‘  (""  " TTTp - )(sec20+l)xr,  + 


2+2 p  _  2+2p 

(X2  "  X1  V.  2  .  (54) 

V  2  +  2p  Jo 


These  four  coefficients,  constant  for  any  given  configura¬ 
tion,  are  specified  by  purely  geometric  conditions,  namely, 
the  cone  angle  0  ,  pivotal  axis  location  xq  ,  the  co¬ 
ordinates  x |  and  x2  which  specify  the  ablation  region 
and  the  exponent  p  which  controls  the  x-dependence  of  the 
ablation  distribution. 


After  car  y  ig  out  the  integrations  in  Eq.  (50)  and 
introducing  the  coefficients  Cg  ,  C^  ,  Cg  and  Cg  there 
is  obtained: 


cm(M)  =  - 


4r| 


"it* 

tan^0 


i  X 


(a  C, 


W 


M— 

tan0 


C5> 


~  „  2  Q  2  i  A 

2t1  Kr  6e 

tan  0 


(a  C4  + 


_ 9^ _  r  ) 

2U  tan  6  V 


(55) 


After  using  the  relations 

•  ^  2  j  ^ 

e1  =  cosA  +  i  sinA  ,  el  -  cos  2A  +  i  sin  2 A  and  q  =  iwa 

In  Eq.  (55)  and  separating  the  real  and  imaginary  parts, 
there  follows 
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C M(p,t)  =  - 


4n  kr0  r 


tan  8 


_a(c3  cosX  “  2lr  tarTe  C5  Sin7^  + 


iauud  ir  .  2U  tan9  r  ,*1 

SU  tan 0  C5  cosX  +  uJB  C3 


2*  Ku  6 


tan  9 


[a(c4  C0S  2X-  2irf^c6 5in  2X1  + 


■JITfanS  (C6  005  2X  +  WTSF5  C4  sin  2X)]  '  !561 


The  pitching  moment  derivatives  are  readily  obtained  from 
Eq.  (56)  by  appropriate  differentiation  giving 


cMa(P,-t) - tan2g  {C3  cosX  ‘  HTTane  C5  Sin>v) 


2t1  K  /e 

- (C,  cos  2 * 

tan^G  * 


T^C6  sin  2*> 


acM(p,t) 


4T1  KR9  f2U»/r  .  u)d  C5cos^\ 

"  "  T77  (¥dKc3Sin*  +  2F~uHT7 


taire 


2ti  KR  9  ,2U 


tan  6 


(|^)(c4sin  2*  + 


uud  C6C0S  2*\ 
7C  tan 9  y 


Eqs.  (57)  and  (58)  express  the  static  and  dynamic  pitching 
moment  derivatives  produced  by  ablation  for  the  general 
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laminar  boundary  layer  case  in  terms  of  the  aerothermo- 
dynamic  quantities  rj  and  ,  the  geometric  coeffi- 

.  ,  and  cone  angle  8  ,  the 
and  the  phase-shift  angle,  X  . 


cients  C, 


'  C4  '  C5  ' 


reduced  frequency. 


uu  d 

7U 


3.4.2  Turbulent  Boundary  Layer 


For  the  turbulent  boundary  layer  case,  m  =  l.o  . 
The  binomial  coefficients  of  Eq.  (48),  when  evaluated, 
give: 


(m+1)  =  2.6  ;  ("J1)  =  .416  ;  l"1*1)  =  -.0112 

2m  =  3.2  ;  f2")  =  1.408  ;  (2|£)  =  -.0113 


By  substitution  of  these  values  in  Eq.  (48)  and  noting 

a  /ae\5 

that,  for  |-^p|  <  1  ,  the  terms  involving  f-g-j  may  be 

neglected  with  respect  to  the  other  terms,  there  is  then 

obtained 

ru/2  1 . 6_  p  i  X  .  a  2n 

I  j ^  =  J  Cp(P)sincp  dcp  =  -5.2titt  K^9  x  e  1+  .12  (-^p )  ) 


,  ,  ,  2p  r,  ,  a  2. 

-  3.2 t|tt  Klt202,2  X  e2'  ae(1  +  -33(-^)  )  (59) 


where  the  subscript  t  is  now  used  to  designate  the  tur¬ 
bulent  boundary  layer  conditions.  By  substitution  from 
Eq.  (59)  into  Eq.  (47),  the  expression  for  the  pitching 
moment  coefficient  produced  by  Elation  becomes: 
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C««P.t>  -- 


5-6  V 

5.2ti  K . .  0  e  p  2  2+t,  2  _1**P_  \  r  ap  2-i  _ 

1  -  J-  (x  sec  9_x  xo/e-  |)  jdx 

*1 


~T 
tan^  0 


2  2-  2  2  i  7\  -  u,n 

3.2ii  K,4  9  e  rx2,  2+2  p  2  _I+Jp.  r  <V  2  _ 

- —V -  J..  l^x  sec  e-x  x Jo  [l+. 331-5)  Jdx 

tan  0  Xj 


(60) 


tq.  (60)  reveals  that.,  unlike  the  laminar  case,  the  moment 

coefficient  for  the  turbulent  case  is  non-linear  in  a 

^  e 

For  small  values  of  -jp  ,  say  |-jpl  <  the  contribution 
of  the  non-linear  terms  will  be  very  small  and  has  been 
neglected  here.  Under  these  conditions,  there  is  obtained, 
after  introducing  the  expression  for  a  of  Eq.  (4)  and 
the  coefficients  [3,^,6^,  Cg  given  earlier  and  car¬ 
rying  out  the  integration  with  respect  to  x  : 


ci|(p't)  =- 


5.2ii  K  1 1 0 1  ’ 6e 1  ^ 

- * - 

tanz  0 


+ 


qd 

TOtanT 


3 . 2t|  Ku202*292'* 
tan  'G 


(“  C4  +  Tirfa KG  Ce) 


(61 ) 


After  introducing  the  complex  trigonometric  exiressions 
for  e1^  and  e2'^  and  using  q  =  iuia  into  Ea.  (61) 
and  separating  the  real  and  imaginary  parts,  there  fol¬ 
lows  : 
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5.2ti  K . ,  0  ‘  d  _  ,  ,  . 

C  ,(M)  =  ^ rQ  ”  [a(c3cosX  -  ^nj^Q  C5sinx)  + 


iuuad  !  ,s  .  2Utan0  ^  ,  M 

TUTin^  ^5C0Sa  +  — 073 —  C3s,n>Vj 


illV—  K 

tan^G  L  v 


C4cos  2X  '  2irfaV'5  C6sin  2V  h 


iuuad  (r _ ,  2Utan0  r .  ; . 

2lft  an  0  lC6  2X  ’  uud  C4°  2>vJ 


The  pitching  r-'ent  derivatives  are  obtained  from  Eq.  (62) 
by  appropriate  differentiation.  Thus 


5 . 2t|  K  0  1  •  ^ 

CMa(P,t)  =  ui?9  (C3cos?'  -  2Utan0  C5sin,>) 


3-21  K , t2  e2 ' 2 

tan^  8 


uud  C6Sin  2X' 


(C4cos 


5. 2rj  K..0 

c  o,  t )  = - ii- 

q  an2  8 


(2lVc  -in'  i  ujd  C5C0SS 

uud  \  3s  n?v  +  7U  TarT^i 


<5  v  2  q2  •  2 

3  .  /LX\  K , ,  8 


It  9  ’  ,2U  Ur  .  .  uud  C6co"  2X' 


tan  8 


/21m /V  „  ,  uud  •'6UUJ 

Hid  \C4S  n  2*  +  7[j  — tinre — /  •  164 


Eqs.  (S3)  and  (64)  express  the  static  and  dynamic  pitch¬ 
ing  moment  derivatives  for  the  general  turbulent  boundary 
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layer  case  with  the  additional  assumption  of  negligibly 
small  nonlinear  terms  in  a  .  They  are  similar,  in  form, 
to  the  corresponding  equations  for  the  laminar  case,  !.e., 
Eqs.  (57)  and  (58).  For  the  same  geometric  configuration, 
the  pitching  derivatives  will  differ  in  magnitude  for  the 
two  cases  on  account  of  the  d  i  f  f  erer.^es  in  the  valuesof 


K 


U 


and  ,  the  constant  multiplying  factors,  the  ex¬ 


ponents  of  0  and  the  coefficients  C, 


and  C 


6* 


3.4.3  Relations  Between  the  Dynamic  and  Static  Abla- 
t i on  Der i vat i ves 


It  is  useful  to  examine  the  relationship  between  the 
dynamic  and  the  static  pitching  derivatives  produced  by 
ablation.  The  values  of  0)/C^  O)  may  be  calculated 

C]  ct 

by  dividing  Eq.  (58)  by  Eq.  (57)  for  the  laminar  case  and 
Eq.  (64)  by  Eq.  (63)  for  the  turbulent  case.  The  result¬ 
ing  expressions  for  the  general  case  are  somewhat  unwieldy. 

K 


However,  i*  ider  rather  realistic  conditions  when  K|  is 
small  (say,  of  order  10  M  and  the  reduced  frequency 

is  low  (say,  of  order  10  )  then  the  terms  !n  K,^  produce 

a  relatively  negligible  contribution  and  tne  relations  may 
be  approximated  as  follows: 


For  the  laminar  case: 


and 


4ti  Ki,9 

Ci.  0,-t)  w  -  - * —  Co  cos^ 

a  t  an2  0  3 


(65) 


4rj  K  i ,  0  9I,  f  mH  CrCosXk 

CM  0,t)  * - -7—  (¥j)(c3sin7v  +  70-fanr)  ‘ 

q  an  0 


(66 ) 
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If,  in  addition,  the  following  inequality  applies 

||u  tan  X|  »  |r-f|nl| 

then,  a  very  simple  relationship  between  the  dynamic  and 
the  static  pitching  derivatives  is  obtained,  namely 

CM  ,PI 

.  (701 

Ma 

Eq.  (70)  illustrates  a  characteristic  effect  of  the  phase- 

shift  angle  X  .  For  a  phase  lag  (A  <  0)  with  0<  |a|  < 

the  effect  of  ablation  on  the  dynamic  derivative  will  be 

of  opposite  sign  to  that  on  the  static  derivative.  It  also 

2U 

shows  that  the  quantity  >  the  indorse  of  the  reduced 
frequency,  acts  as  an  "amplification  factor"  for  the  dy¬ 
namic  derivative  with  respect  to  the  static  one. 
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3.4.4  Relations  Between  the  Ablating  and  Nonablating 
Pi  t  c  h  i  n  g  Moment  Derivatives 

it  is  of  interest  to  relate,  parametrically,  the 
static  and  dynamic  pitching  moment  derivatives  contri¬ 
buted  by  ablation  to  the  corresponding  values  without 
ablation.  Attention  will  be  restricted  to  the  special 
cases  given  by  Eqs.  (65  -  68}.  The  small  contributions 

of  the  spherical  nose  to  C..  (o)  and  C,,  (o)  will  not 

Ma  Mq 

be  considered.  From  Eqs.  (65)  and  (13)  one  obtains,  for 
the  laminar  case: 

C|Yi  O,^)  q 

■C°  |o)  =  2  K,,cos  ’'TT  •  l71) 

‘•Mam  ]t  L1 

For  the  turbulent  case,  Eqs.  (67)  and  (13)  yield 

CM  lfJ’t)  6  C 

— C-M  (or  -  2-6  Klte  cos  •  1721 

Equations  (71)  and  (72)  illustrate  the  relative  importance 
of  the  different  factors  contributing  to  the  fractional 
static  pitching  moment  derivatives,  namely,  the  quantity 
,  the  geometric  coefficients  and  ,  the  cone 

angle  0  and  cos X  .  For  small  phase-shifts,  of  the 
order  of  +10  degrees,  cos>>  «  1  so  that  the  fractional 
static  moment  derivative  is  insensitive  to  phase-shifts 
us  regards  both  magnitude  and  direction.  Further,  since 
^3 

t—  may  be  positive  or  negative  depending  on  the  ablation 
^  1 
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and  cone  geometries,  the  ablation  contribution  to  C.. 

Mqi 

namely,  C,,  (P)  may  be  stabilizing  or  destabilizing. 
wa 

Turning  our  attention  to  the  fractional  dynamic  pitch¬ 
ing  derivatives,  then  Eqs.  (66)  and  (14)  give 

cM  (P,t) 

T?ToT  =  2  KHtan0(ii)(sinA 
Mq 

Eqs.  (68)  and  (14)  yield  for  the  turbulent  case: 

CM 

T^ToT  “  2~6  Klte'6tane,l3,(si"x 


C3 


uud 

2U 


cosA 

t  a  ft  C 


ir2) 


(74) 


Eqs.  (73)  and  (74)  show  that  the  fractional  dynamic  pitch¬ 
ing  moment  derivative  depends,  as  before,  on  the  quantity 

2D 

K|  and  also  on  0  ,  s i nX  ,  cosA  ,  — ^  ,  and  the  ratios 

C3  S 

-x—  and 

l2 

If  the  quantities  involving  -jg  in  Eqs.  (73)  and 

(74)  are  negligibly  small,  relative  to  the  other  term 

2U 

and  if  we  note  that  ,  6  ,  ,  and  C2  are  nonnega¬ 

tive,  it  becomes  apparent  that  tue  dynamic  derivative 
due  to  ablation,  Cu  (P)  ,  will  be  of  the  same  polarity 

Mq 

as  C..  (0)  ,  i.e.,  damping  whenever  C^sinT^  is  positive. 

Mq  J 

Therefore,  positive  values  of  coupled  with  phase 

lead  or  negative  values  of  coupled  with  phase  lag  will 
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produce  the  same  result,  namely,  an  increased  damping. 

On  the  other  hand,  when  and  s ; nA  are  of  opposite 

sign,  the  dynamic  effect  of  ablation  will  be  uestabiliz- 
i  ng. 

3.5  Effect  of  Ablation  on  the  Local  Friction  Coefficient 


The  gases  generated  as  a  result  of  ablation  and  in¬ 
jected  into  the  boundary  layer  may  produce  a  significantly 
large  change  of  the  local  skin  friction  coefficient.  This, 
in  turn,  produces  a  corresponding  change  in  the  pitching 
moment  coefficient  and  the  pitching  moment  derivatives, 

designated  here  by  C^(f)  >  (f)  ,  and  (f)  . 

a  q 

The  dependence  of  the  local  skin  friction  coefficient 
on  ablation  or  blowing  will  involve  the  local  mass  transfer 
parameter  as  well  as  other  flow  characteristics  such  as  the 
wall  to  free  stream  temperature  ratio  and  the  mol  cular 
weight  ratio  of  the  free  stream  and  vaporized  gases. 

Li  and  Gross^  have  proposed  the  following  expres¬ 
sion  for  a  compressible,  laminar  boundary  layer  with  blow¬ 
ing  or  suction 


(75) 


where  and  are  the  skin  friction  coefficients 

0  with  and  without  blowing,  respectively; 
A  :s  a  positive  quantity,  fixed  for  a  given  fluid; 
T  and  Tq  denote  the  uniform  w_  1 1  temperature  arid 
the  free  stream  stagnation  temperature, 
respect i ve ly; 
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and  K  is  a  nond imens i ona  1  mass  transfer  parameter  such 
that  K  <  0  fcr  blowing  and  K  >  0  for  suction, 

12 

Gross,  Hartnett,  Masson  and  Gazley  have  examined 
the  available  analyses  of  the  binary,  laminar  boundary 
layer  and  concluded  that,  for  engineering  design  calcula¬ 
tions,  the  skin  friction,  in  the  presence  of  mass  transfer 
cooling  with  a  foreign  gas,  may  be  approximated  quite  well 
by  the  expression 


(765 


where 


and 


m 

— -  represents  the  ratio  of  the  molecular  weights 
m 

c  for  air  and  the  coolant  gas; 

F  is  the  nond imens i ona 1  mass  transfer  parameter; 
Rex  is  the  local  Reynolds  number 
C*  is  the  Chapman-Rubes i n  parameter,  calculated 
at  the  reference  temperature. 


For  the  turbulent  boundary  layer  case,  Danberg, 

1 8 

Winkler  and  Chang  obtained  experimental  results  for 

the  effects  of  air  transpiration  and  wall  temperature 

on  skin  friction  on  a  flat  plate  at  a  Mach  number  of 

6.7  indicating  a  strong  dependence  of  the  skin  friction 

coefficient  on  mass  transfer  even  for  the  relatively 

low  mass  transfer  coefficients,  F„  =  ,0025  Testsen- 

e 

compassing  the  following  range  of  parameters 


0  £  F  *  .002  5 
e 


0  * 


2F 

e 


<  4 


.2 
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yielded  the  following  simple  expression: 


=  1  -  .22 


Equations  (75),  (76),  and  (77)  are  all  of  the  generalized 
type: 


=  1  -  hF 


where  the  nond imens i ona 1  factor  k  is  determined  by  the 
relevant  parameters  of  the  problem  at  hand.  It  replaces 
the  quantity  ./"^wN  cf  Eq.  (75),  the  more  general  quan- 
tity  'V 


ty  cf  Eq. 


of  E  (76)  or  the  quantity  of  Eq.  (77)  for  the 


turbulent  boundary  layer  case,  o 

If  one  refers  to  Figure  1  and  focuses  attention  on 
the  elemental  normal  and  axial  forces  due  to  friction, 
dN(f)  and  dX(f)  ,  respectively,  acting  on  an  element 
of  surface,  dS  ,  on  the  cone,  one  may  write 

a N  (  f )  =  qe  Cf  si.nG  s  i  ncp  dS 


wnere 


dX(f)  =  q  C,  cosG  dS 
e  f 

i  q.  -  '/2  p„u  2 
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3 . 6  Pitching  Moment  Derivatives  Due  to  Changes  in  Skin 
Friction 


The  elemental  pitching  moment,  dM(f)  ,  about  the 
pivotal  axis  at  xq  can  be  expressed  by 


dM ( f  )  =  dX{f)  r  s  i  ncp  -  dN(f)  (x-x  ) 


179) 


Noting  that  r  =  x  tan9  and  dS  =  dlPdx 

c0’jW'  cos^a 


then,  by  substituting  for  d  X { f )  ,  dN(f)  ,  r  and  dS  in 
Eq.  (75)  and  simplifying,  there  is  obtained 

2 

dM { f  )  =  q  tan  0  Cfxx  s  i  ncp  dcpdx 
0  r  0 

The  correspond i ng  pitching  moment  due  to  friction,  M(f)  , 
is  obtained  by  integration  over  the  cone  surface 

rL  »tt/2  ->  pL  pir/2 

M{f)  =  2J  J  dM { f )  dcpdx  =  2q  tan  8J  j  C^sincp  xxQ  dcpdx 
Lj  -tt/2  L 1  -tt/2 

(80) 


V/hen  we  introduce,  as  before,  the  nond  i  mens  i  ona  1  quanti¬ 
ties  x  ,  x  ,  L,  ,  use  the  cone  base  area  AD  and  base 
diameter  d  as  reference  area  and  moment  arm,  respectively, 

and  refer  the  moment  coefficient  to  the  free  stream  dy- 

2 

namic  pressure,  qm  =  l/2  ,  then  the  moment  coeffi¬ 

cient  due  to  friction  will  be  given  by 


M  ( f  ) 


pTt/2 


Me  r  r  ' 

rffan^  JT  J_TT/2Cf  xxo  sincP  dcpdx 


(81  ) 
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By  virtue  of  E  q  (81),  the  pitching  moment  coefficient  due 
to  friction  can  be  calculated  once  the  functional  relation 
between  the  local  skin  friction  coefficient,  ,  and  the 
variables  cp  and  x  is  known.  Equation  (78)  peonits  the 
contribution  to  the  local  skin  friction  due  to  ablation, 
designated  here  by  AC^  to  be  written  as  follows: 


AC 


f 


k  F  C, 
e  f 


(82) 


The  corresponding  pitching  moment  coefficient,  ( f ) 
may  be  obtained  by  means  of  Eq.  (81)  but  with  AC^  in¬ 
stead  of  and  with  the  limits  of  integration  of  x 

determined  by  the  ablation  region.  Hence 


CM(f) 


q  ,*2^/2 

q-T,tane  *° 


s  i  ncp  dcp  dx 


(83) 


In  Eq.  (83),  k  and  Cf  are  independent  of  the  angle 

'  o 

cp  so  that  one  may  write 


CM(f> 

qe 

,x2 

_  [.tt/2 

Cr  x  x  F  s  i  ncp  dcp  dx 

fo  0  -tt/2  e 

(84) 

q  ntanS  J 

K 

xi 

if  we  use 

for  F 

e 

the 

expression  obtained  from  Eqs. 

(36) 

and  (40), 

namely 

then.  Eq.  (84)  becomes 

q  K.gV*  rx2  1+p  rn/2  m 

CM(f }  =  ’  Tq^TanT*  J-  HCf  x  xoJ  /9y  'sincp  dcp  dX  *  (85) 

x  j  o  -tt/2 
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a  s  i  ncp 


We  restrict  consideration  to  the  case  |-£— ~| 
use  the  binomial  series  expansion 


m  /  aes  i  ncpxm  .a  si  ncp,  n 

y  =  l1  -  -h—j  -X  <-'>  O  (-V-) 


Considering  first  the  last  integral  in  Eq.  (85): 


'2  (-,,n  <n>  (£)"  l^‘?>n+'  d* 


It  Will  vanish  for  zero  and  even  values  of  n  while 
f  or  odd  va 1 ues,  it  y  i  e  1  ds 


=  -  T  -  _£  \  3  _  /ae\3  ,m »  5  /ae\° 

2  '  2  9  *  {3}  e^VTJ  ‘  l5}  T6n(~)  " 


,a  5 


By  subst  i  tut  ion  into  Eq.  (35)  there  is  obtained  the  follow¬ 
ing  genera!  expression  for  the  pitching  moment  contribu¬ 
tion  due  to  friction 


Vfi  * 


qe  K|9Wda  i-x2 


a  t  a  n 0 

GO 


r  *  i+p  _  a 

I.  *  c,  x  ;  r.  £  « 

Jz  f  oL  7  0 


<"> I  ffl  - 


a  ,3 


<5>  f-TS1  (^)5j  95 


3o6.1  Laminar  Case 

For  the  laminar  boundary  layer  case,  m  -  1  and  only 
the  first  term  inside  the  square  brackets  remains  so  that 
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1/  i  A  x~  ,  . 

qpKi^e  r  2  -,+p- 

~  fq  atanS  J-  HiCf  x  x0aedx 

®  X  1  0  J  -v 


If  one  uses  for  the  very  genera)  expression,  derived 

from  Eqs.  (76)  and  (78),  namely 

m  1/3  /TTT 
,m,\  J  x 

K.  =  2.08  (~r)  — -  (81 

-c  \m  '  /7rr 


and,  for  C. 


the  general  expression  given  in  Refer¬ 


ence  44  modified  to  apply  to  the  case  of  the  cone,  i.e.. 


Cf  -  *  * -£§44- (ft) 

0,1  /  x  e 


where  the  subscripts  *  and  e  denote  the  reference 
temperature  and  gas  temperature  at  the  edge  of  the 
boundary  layer,  respectively,  and  e  denotes  the  ex¬ 
ponent  in  viscosity  power- lew  relationships.  Then,  by 
combining  Eqs.  (88)  and  (89),  it  is  apparent  that  the 

product  will  not  involve  the  local  Reynolds' 

o ,  "t 

number  and  may  therefore  be  taken  outside  the  integral 

sign.  Therefore,  after  using  Eq.  (4)  for  ag  and  the 

abbreviation,  -  h^C r  ,  Eq.  (87)  becomes: 

c  1  o,l 


.  qo  e1  *  x2_  i+p_ 

Cjf,t)  =  -2  -  X  X 

to  qm  2atan8  J-  o 

•oo  y 


rX2  i+p  r  qd(x-x  )-,  _ 

J-  x  xoLa  +  mzzs  -jdx 
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2+p  2+p 

Let 

‘7  '  ft 

(91 ) 

and 

_  3t p  3+p 

2+p  2+p 

o  _  fX 2  '  X1  )z  _ 

(*2  '  X1  )r.  2 

(92) 

8  V  T^Tp /  o 

These  two  coefficients  are  determined  by  purely  geometric 
conditions  and  are  therefore  constant  for  a  given  configura¬ 
tion. 


After  integration  of  Eq.  (90)  and  the  introduction 
of  the  coefficients  Cy  and  r.g  ,  the^e  is  obtained 


CM(f,4) 


KRK2i 

Qco  2 a an 0 


+ _ 

+  2Utan0 


(93) 


After  using  the  relations  e1^  =  cosA  +  i  s >  nA  and 
q  =  iuua  in  Eq.  (93)  and  separating  the  real  and  imagi¬ 
nary  parts,  the  following  expression  is  obtained  for  the 
moment  coefficient: 


CM(f,4) 


%  KuK2t  r 

q_  2  at  an  9 


[°(C7  C 


cosA  - 


C  s 

2Utan9 


i  nA) 


2U  (“TiTTS —  +  I?  c7  sin70] 


(94) 


The  pitching  moment  derivatives  due  to  the  frictional 
effects  of  ablation  readily  follow  from  Eq.  (9^)  by  ap¬ 
propriate  differentiation 
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c..  lf,D 


%  *UK2 1 

2  a  t  a  n  8 


cosX  - 


md  C8sin\ 
21)  t  a  n  §  J 


and 


Cu  {f’l) 

Mq 


%  KHK2 l 

qm  2atan6 
00 


s  i  nX 


,  CqCOsX-, 
.  u)d  8  I 

TU  tan0  J 


(95) 


(96) 


3.6.2  Turbulent  Case 

For  the  turbulent  boundary  layer  case,  m  =  1.6  and 
the  binomial  coefficients  in  Eq.  (86)  give 


m  *  1.6  J  (3)  =  -.064  ;  and  (™)  =  -.0108 


By  substitution  of  the  i  va lues  into Eq.  (86)  there  is 
obta  i  red 

1/  q  •  6  I  X  X/n  1  , 

q_  ^ I f 8  6  p  2  _  '  +  P_  r 

=  ’  q  “Fan" 0~  J-  K2tX  xo["‘8  ae  + 

»  X1 


x2 


a  v4- 


•024(t)  +  -003  (f)  Y* 


(97) 


Unlike  the  laminar  case,  the  moment  coefficient  produced 

by  turbulent  friction  is  not  linear  in  ct  .  However, 
a  e 

for  i~(  <  1  ,  the  nonlinear  terms  are  quite  small  so 

w 

that  only  the  first  term  inside  the  square  bracket  of 
Eq.  (97)  need  be  retained.  Noting  also  that  the  quantity 
l<2j.  is  independent  of  x  ,  then  after  s  imp  1  i  f  i  cat '  on  and 
integration  of  Eq.  (97)  and  the  introduction  of  the  coef¬ 
ficients  C 7  and  Cg  ,  there  is  obtained: 
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qp  KitK?te 


.6 


e'X[a  C7  +  2Utan0  C8J 


(98) 


Following  the  procedure  described  earlier  to  resolve  Eq. 
(94)  into  its  real  and  imaginary  pjrts  yields  the  follow¬ 
ing  expression  for  the  turbulent,  frictional  pitching 
moment  coefficients: 


CM(f^ 


,  K , ,  K„  ,  0 
.  8  *  2t 


Q» 


.6 


t  an  6 


w 


...  j  C q s  i 

C7  cos*  '  TU  "TanT" )  + 


Cgcosx  2U  r  c. 

+  — — r  C-j  s  i  n^  J  J 


tanQ 


it'd 


(99) 


The  corresponding  pitching  moment  derivatives  will  be 

qa  K.-K^.Q*6  ,  Cnsin\ 

CMa(fjt)  =  *8  dl  "TTanG  (C7  cos *  '  TIT  “ TTJ 

and 


CM  (f,t)  = 


■B  qeKltK?t9 
a  qm  tan0 


.  6 


,aU)[c7  . 


(  100) 


(101) 


3.6.3  Relation  Between  the  Dynamic  and  Static  Fric¬ 
tional  Pitching  Moment  Derivatives 

The  relation  between  the  dynamic  and  the  static 
derivatives  contributed  by  frictional  effects  may  be 
calculated  from  Eos.  (96)  and  (95)  for  the  laminar  case 
and  Eqs.  (101)  and  (100)  for  the  turbulent  case.  In 
general,  the  following  inequality  holds: 

^UtariQ  CS  sin?sl  <<:  lC7  cos^l 
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Therefore,  Eqs.  (96)  and  (95)  and  also  Eqs.  (101)  and 
(100)  yield 

Cu  (f)  c 

Mq  _  2U  .  .  .  C8 

CM  (f )  “  wd  tanX  ITJTane 
Ma 

Eq.  (102)  is  valid  for  both  the  laminar  and  turbulent 

cases;  it  is  similar,  in  form,  to  Eq.  (69)  but  with 
co  C5 

replacing  7 r- 
U7  L3 

With  the  further  restriction  that 

liH  tan>'l  >:>  iTTfanSl 

Then,  the  very  simple  relationship  between  the  dynamic 
and  the  static  pitching  derive  ives  found  to  apply  to 
the  inviscid  effect  is  also  valid,  namely, 

CM  (f! 

q  2U 

Cu  If)  “Sd  tanl  ■ 


(102) 


(103) 


-45- 


FluiDyne  Engineering  Corporation 


4.  AVERAGE  HEAT  AND  MASS  TRANSFER  RATES 


In  Section  3,  explicit  expressions  were  given  for 
the  local  heat  transfer  on  the  cone  surface,  t)(c)  ,  for 
both  the  laminar  and  turbulent  boundary  layer  conditions. 
Through  the  concept  of  the  effective  heat  of  ablation. 


the  corresponding  local  rates  of  mass  ablation. 


m  are  readily  calculated: 
w 


In  order  to  calculate  the  overall  rates  of  heat  transfer 

•  • 

and  mass  ablation,  Q(c)  and  Mw  ,  respectively,  it  is 
only  necessary  to  integrate  the  local  values  over  the 
ablating  surface  on  the  cone.  Therefore 


0(c)  =  J  Mc)dS 
JS 

and 


rx2pn/2 

2 

x ,  -tt/2 


)x  sin0  dtp  dx 

- rj—— — i - 

cos*"  0 


2  f  2  r"/2  fr(c)x  si  n9  dtp  dx 
^Jx1  -tt/2  Heff  cos2  8 


(104) 


( 105) 


By  substituting  the  general  expression  for  Cj  ( c ) 
given  by  £q.  (29)  and  introducing  nond i mens i ona 1  quanti¬ 
ties  as  before,  there  is  obtained: 


.  n.  2+P„  •  a  [.tt/2  1+p 

0(c)  =  2i - CSLne  r  '  X 


(106) 


a  s  i  ncp 

Noting  that  =  0mym  where,  as  before  y  =  1 - ^ - 

then,  using  the  binomial  series  expansion,  there  is  ob- 
ta i ned  f r om  Eq.  { 104  ) : 
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Q(c)  = 


„  2+pr  Qrn  .  0  rx2  pTr/2  “  n  m  1+P  aa  n 

2L - Cl  sm8  f  J  (-|)"(")5  (-g^)  (si  ncp)n- d<pdx 

cos  0  -n/2 


n=0 


(107) 


The  integrals  involving  odd  powers  of  sincp  will 
vanish  for  the  present  limits  of  cp  and  the  expansion 
of  Eq.  ( 107)  yields: 

Q(C)  =  liZ^.sin.e  J'2;Hp[t,  +  (”)(^)2f+  C4)&>4jgn  ■■■ 


CCS  0 


4. 1  Laminar  Boundary  Layer 


(108) 


For  the  iaminar  case,  m  =  1  and  Eq.  (108)  is  readily 
integrated  and  yields 

2+P, 

'  . .  p  t-  _  ,TK 

J  -  *  tt  dx 


Q^c) 


2LtTfJC^0  sine  ‘'2  _  1+p 


T 

cos4  0 


2tt  L2+pC^0  siri0  ~  xi  ^ 

- ^ -  - 

cos  0 


2+p  _  2+p 

/x2 

l - F+~p - J 


2tt  L2+pC^0  sin0  C? 

- - -  — 

cos  0  x„ 


(109) 


The  overall  rate  of  mass  ablation  will  be  given  by  the 
express i on 
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rjt(c)  2tt  L^+pC^6  sine  C? 


:eff 


H  ,,  cos  0  x„ 
eft  o 


(no) 


4.2  Turbulent  Boundary  Layer 


For  the  turbulent  boundary  layer  case,  m 
values  of  the  binomial  coefficients  are: 


=  1.6  .  The 


(*)  =  .48  ;  (J)  =  .0224  . 


Eq.  ( 108)  then  yields: 

o-rr  i  2+P r  n  ^  : 


2tt  Lz'rMC+e' ’“sinQ  r  2  l+pr  a  2  a  4 

0 1  ( c  )  = - yi -  J.  x  [l  +  .24(-g£)  +  .  007  (-jp) 


cos  8 


(111) 


,a  .  4 


Since  j-jpj  <  1  ,  the  term  involving  is  negligibly 

small  and  will  be  deleted.  Eq.  (Ill)  shows  that  the  over¬ 
all  rate  of  heat  transfer  for  the  turbulent  boundary  layer 
case  involves  some  dependence  on  the  angle  unlike  the 

laminar  case  as  given  by  Eq.  (109)  which  exhibits  no  such 
dependence. 

If  the  following  inequality  holds: 


,a  x9 

| .  2  4  \jj~j  !  <<  1  ,  t  h  e  n 


Q  t  ( c )  m 


2n  L2+pCt01,6sin0  Cy 


cos  8 


(112) 
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and  the  corresponding  rate  of  mass  ablation  is  given  by 


M 


w,t 


2tt  L2+pCte1,6sin9  C? 


H 


eff 


T~ 

cos  0 


(112) 
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5.  ILLUSTRATIVE  EXAMPLE 

In  order  to  illustrate  the  application  of  the  fore¬ 
going  analysis,  a  hypothetical  body  was  chosen  and  the 
flight  configuration  selected  was  such  as  to  reproduce 
the  conditions  typical  of  laminar  and  turbulent  boundary 
layers  successively  during  the  re-entry  phase. 

The  heating  rate  experienced  during  re-entry  depends 
primarily  on  the  time  history  of  velocity,  altitude,  and 
attitude.  For  the  purpose  of  this  example,  attention  is 
focused  primarily  on  the  changes  in  the  dynamic  pitching 
derivatives  due  to  ablation  and  a  simplified  approach, 
assuming  constant  missile  parameters,  is  adopted.  The 
flight  conditions  are  based  on  some  assumed  but  repre¬ 
sentative  altitude,  velocity,  and  oscillation  conditions 
rather  than  on  a  calculated  trajectory  involving  continu¬ 
ous  changes  in  altitude,  velocity  and  atmospheric  condi- 
t i ons. 

The  geometric  characteristics  of  the  hypothetical 
re-entry  body  are  as  follows: 

Length  L  =  10  feet 
Cone  Half-Angle  9  =  10° 

Cone  Base  Radius  r^  =  1.745  feet. 

Bluntness  ratios  of  5=0  (pointed  cone)  and  .2  are 
considered.  The  center  of  gravity  location  xQ  is 
assumed  to  range  between  .3  and  .8  with  .1  intervals. 

A  phase-shift  angle  of  X  =  -.349  rad.  (20°  lag)  is 
assumed . 
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5. 1  Laminar  Boundary  Layer  Case 

The  typical  flying  conditions  chosen  for  this  case 

are : 


Altitude  h  =  150,000  feet 
Velocity  U_  =  20,000  ft/sec 

Free  Stream  Density  pa  =  3.05  x  10”  slug/cu.  ft. 
Oscillation  Frequency  f  =  3  cps 
Initial  Angle  of  Oscillation  aQ  =  5° 

Velocity  of  Sound  aQ  =  971  ft/sec 
Ratio  of  Specific  Heats  y  =  1.35  . 

Calculations  are  made  for  three  arbitrary  ablation  con¬ 
figurations,  namely,  full  length  ablation  ( x  ^  , 

x2  =  1),  ablation  from  the  front  half  only  { x  ^  =  t1  , 

'<2  =  .5)  and  ablation  from  the  rear  half  only  (x^  =  .5  , 

X2  =  ^ ) • 

The  value  of  the  pressure  correction  factor  "n  for 
the  present  case  as  given  by  Eq.  (2)  is  1.057.  The  re¬ 
sults  of  the  calculations  for  the  laminar  boundary  layer 
case  are  illustrated  in  the  following  figures. 

The  geometric  ablation  coefficients,  ,  C^  ,  Cg  , 
and  Cg  ,  calculated  for  the  three  ablation  configurations 
and  different  pivotal  axis  location  are  illustrated  in 
Figures  6,  7,  6,  and  9,  respectively.  It  may  be  noted 
that  C0  and  C.  which  are  linear  in  £  may  be  posi- 
t i ve  or  negative  whereas  the  parabolic  coefficients  Cc 

v 

and  Cg  are  always  positive. 
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The  geometric  coefficients  associated  with  the  fric¬ 
tional  effects  are  illustrated  in  Figures  10  and  11.  Here, 
the  linear  coefficient  C -j  is  always  positive  while  the 
parabolic  coefficient  Cg  may  be  positive  or  negative. 

The  static  and  dynamic  derivatives  due  to  ablation, 

(0)  and  (0)  are  calculated  from  Eqs.  (57)  and 
q 

1,  respectively,  for  the  three  ablation  configurations 
;  -he  results  are  illustrated  in  Figures  12  and  13.  The 
calculations  are  made  for  a  phenolic  nylon  ablator  with 
th^  following  rather  arbitrary  assumptions: 

H  ff  =  8000  BTU/lb 
a  =  10 


The  calculated  value  for  for  this  case  was  found 

to  be:  K ^  . 0157. 

The  static  and  dynamic  derivatives  due  to  friction, 

(f)  and  (f)  calculated  from  Eqs.  (93)  and  (94), 
ct  q 

respectively,  are  illustrated  in  Figures  14  and  15.  For 
these  calculations,  an  arbitrary  constant  value  of  .44 
was  used  for  the  quantity  K ^ ^  >  it  is  based  on  Refer¬ 

ence  18. 


The  resultant  static  and  dynamic  pitching  moment 
derivatives  for  the  ablating  cone  are  calculated  by  sum¬ 
mation  of  the  various  contributions  and  the  results  of 
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the  summation  are  illustrated  in  Figure  16  for  the  static 

derivative,  Cu  and  in  Figure  17  for  the  dynamic  deriva- 
m 

tive,  Cu  . 

Mq 

The  relationship  between  the  ablating  and  nonablating 
dynamic  derivatives,  as  given  by  the  simplified  Eq.  (73), 
is  illustrated  in  Figure  18.  It  is  noteworthy  that,  for 
typical  location  of  the  pivotal  axis,  say,  xQ  =  .65  , 
full  length  ablation  on  the  pointed  cone  contributes  a 
dynamic  pitching  derivative,  C..  (0)  of  the  same  sign  and 

Mq 

about  the  same  amplitude  as  the  aerodynamic  damping  deriva¬ 
tive.  Forward  and  aft  ablation  will  contribute  dynamic 
derivatives  of  magnitudes  two  to  three  times  the  nonabla¬ 
tion  value  and  such  that  forward  ablation  increases  the 
dynamic  damping  while  aft  ablation  destabilizes  the  oscil¬ 
lation. 


The  relative  effect  of  the  contribution  of  friction, 
as  modified  by  ablation,  to  the  dynamic  pitch  damping 
derivative,  i.e.,  the  ratio  [f)/C^  (o)  was  calculated 

q  q 

from  Eq.  (102)  and  its  approximate  form,  Eq.  (103),  and  the 
results  are  illustrated  in  Figure  19.  It  is  apparent  that 
these  contributions  have  the  same  sign  as  the  nonablation 
dynamic  derivative,  i.e.,  they  are  always  damping.  For 
the  present  example  and  typical  pivotal  axis  location,  the 
contributions  range  between  10  and  30  percent  for  the  dif¬ 
ferent  ablation  configurations. 


It  is  instructive  to  examine,  in  turn,  the  effects 
of  changes  in  the  frequency  of  oscillation,  f  >  or  in  the 
amplitude  of  the  phase  lag,  ?v  ,  on  the  resultant  dynamic 


derivative,  C 


M. 


Figure  20  illustrates  the  results  of 
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calculations  with  the  three  chosen  ablation  configura¬ 
tions  and  the  same  phase  lag  of  20°  but  with  assumed 
oscillation  frequencies  of  2,  3,  5,  and  10  cycles  per 
second.  The  change  of  frequency  alters  the  slope  of 

the  lines  for  Cu  versus  x  and  the  lower  the  fre- 

o 

Q 

quency,  the  steeper  the  lines. 

Figure  21  shows  the  effects  of  changing  the  phase 
lag  angle  while  maintaining  all  the  other  factors  con¬ 
stant.  Here,again,  the  effect  of  increasing  the  phase 
lag  is  to  steepen  the  slope  of  the  lines  of  versus 


Figures  20  and  21  show  a  striking  similarity  in  the 
variation  of  the  resultant  dynamic  derivative  for  the 
case  of  a  fixed  phase  lag,  (X  =  -20°)  with  different  fre¬ 
quencies  (f  =  2  to  10  cps)  and  that  of  a  fixed  frequency 
(f  =  3  cps)  and  varying  phase  lag  (X  =  -5  to  -40  degrees). 
This  is  due  to  the  fact  that,  for  the  present  example, 
the  contribution  of  ablation  to  the  dynamic  derivative 
is  essentially  i ..  direct  proportion  to  ( s  i  n  >  )/ f  ,  all 
other  factors  remaining  unchanged. 

5 . 2  Turbulent  Bou.  dary  Layer  Case 

The  same  geometric  ablation  configurations  were 
used  and  the  representative  flying  conditions  selected 
for  this  case  were: 

Altitude  h  =  100,000  ft 

Velocity  =  17,500  ft/sec 

Free  Stream  Density  p9  =  3.31  x  10"5  slug/cu.  ft 
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Oscillation  Frequency  f  =  3  cps 
Initial  Angle  of  Oscillation  aQ  =  5° 
Velocity  of  Sound  aQ  =  971  ft/sec 
Ratio  of  Specific  Heats  y  =  1.35. 


The  same  values  for  H 


ef  f 


a  ,  h  /h„  and  were 


w 


used  to  afford  a  simple  comparison  even  though  these 
values  may  not  be  equally  realistic  for  both  laminar 
and  turbulent  boundary  layers. 


The  value  of  the  pressure  correction  factor,  ti  , 
calculated  for  this  case  was  1.063. 


The  geometric  coefficients  Cg  through  Cg  will 

be  numerically  a  little  different  from  the  corresponding 

values  for  lhe  laminar  case  because  of  the  different 

value  of  the  exponent  p  being  -.2  instead  of  -.5  . 

The  variation  of  these  coeffici  nts  with  5  will  have 

o 

the  same  general  characteristics  as  their  laminar  counter¬ 
parts  and  will  not  be  given  in  detail.  The  static  and 

dynamic  derivatives  due  to  ablation,  Cu  (fJ)  and  Cu  (0) 

cx  q 

calculated  from  Eqs.  (63)  and  (64),  respectively,  are 
illustrated  in  Figures  22  and  23.  The  calculated  value 
for  was  found  to  be  .0450.  The  dynamic  derivative 

due  to  friction,  Cu  (f)  calculated  from  Eq.  (101),  is 

H 

illustrated  in  Figure  24. 

The  resultant  dynamic  pitching  derivative  is  illus¬ 
trated  in  Figure  25,  while  Figures  26  and  27  present  the 
variation  of  the  ratios  of  (0)/C..  (o)  and  C^;  (H/C^  (o) 

Mq  Mq  q  q 

with  xQ  for  all  three  ablation  configurations.  These  figures 


-55- 


FluiDyns  Eng'niehinu  Corpopatioh 


exhibit  the  same  general  characteristics  as  for  the  lami¬ 
nar  case.  There  is,  however,  a  significant  difference 
for  the  case  of  full  length  ablation.  The  contribution 
of  ablation  to  the  dynamic  derivative  tends  to  destabilize 

the  oscillation  for  values  of  x  below  .66  for  the 

o 

pointed  case  compared  to  .6  for  the  laminar  case.  The 
practical  significance  of  this  shift  is  discussed  in  the 
following  paragraph. 

Figures  28  and  29  provide  a  graphic  comparison  of 
the  resultant  static  and  dynamic  derivatives,  respectively, 
calculated  for  the  three  different  ablation  configurations 
on  the  pointed  cone  for  the  two  flight  cases  which  illus¬ 
trate  the  laminar  and  turbulent  boundary  layer  conditions. 
From  Figure  28,  it  is  apparent  that  at  any  specified  value 
of  xQ  ,  the  differences  in  the  value  of  the  resultant 
static  derivatives  for  the  laminar  and  turbulent  flow  con¬ 
ditions  or  for  different  ablation  configurations  are  quite 
small.  Figure  29,  on  the  other  hand,  illustrates  very 
significant  differences  in  the  relationships  between  the 

resultant  dynamic  derivatives  C..  and  x  for  the  lami- 

Mq  u 

nar  and  turbulent  boundary  layer  cases.  Specifically, 
with  full  length  ablation,  the  pointed  cone  is  dynamically 
stable  for  x„  >  .56  for  the  laminar  case  and  x„  >  .61 
for  the  turbulent  case.  In  this  practically  important 
range  of  center  of  gravity  location,  the  resultant  dy¬ 
namic  pitching  moment  derivatives  will  have  different 
signs  depending  on  the  type  of  boundary  layer.  This 
feature  may  have  far-reaching  consequences  for  conical 
vehicles  re-entering  the  atmosphere.  At  high  altitudes 
of  the  order  of  150,000  feet  and  flying  speeds  of  20,000 
ft/sec,  the  ambient  conditions  will  be  such  that  the 


-56 


FlujDyns  Engineering  Corporation 


unit  Reynolds  number.  Re/ft  ,  .Ml  be  shout  1.5  x  105 
The  boundery  layer  is  laminar  and  the  resultant  dynamic 
der, vat, ve  will  remain  negative,  i.e.,  dynamically  Cable. 
;  a'tltuiJes  of  about  100,000  feet  and  air  speeds 

°  Z*Vk  '7'500  ft/se'*  th«  “"it  Reynolds  number  is  abet 
*  so  that  the  boundary  layer  becomes  turbulent 
except  near  the  nose,  and  the  dynamic  derivative  becomes 
Pos.t.ve  so  that  the  vehicle's  oscillations  will  diverge. 
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6.  SUMMARY 

A  simple  analytical  method  has  been  presented  which 
permits  the  calculation  of  the  hypersonic  static  ar.d  dy¬ 
namic  pitching  moment  derivatives  for  slender,  pointed, 
or  slightly  bunted  cones  performing  s  i  .»g  1  e-degree-of- 
freedom,  snial  1-angle  oscillations  about  the  center  of 
gravity,  taking  into  account  the  effects  of  ablation 
over  arbitrary  portions  of  the  cone  surface-  The  analy¬ 
sis  is  based  on  Newton’s  Impact  Theory  for  hypei  sonic 
flow  and  is  restricted  to  considerations  ot  rigid  body 
motions  at  low  reduced  frequencies.  Quasi-steady  condi¬ 
tions  are  postulated. 

The  resultant  static  and  dynamic  derivatives  are 
obtained  by  adding  to  the  nonablating  derivatives,  sep¬ 
arate  add'tional  contributions  due  to  local  mass  trans¬ 
fer  effects  and  to  changes  in  the  local  skin  friction 
resulting  from  the  ablation  process. 

The  local  heat  transfer  rates  and  skin  friction  co¬ 
efficients  are  obtained  from  classical  flat  plate  results 
appropriately  modified  to  account  for  the  local  effective 
angle  of  attack  and  the  conical  flow  conditions.  The  con¬ 
cept  of  an  "effective  heat  of  ablation"  is  applied  to  re¬ 
late  the  local  heat  and  mass  transfer  rates. 

The  effects  of  the  ablation  geometry  are  conveni¬ 
ently  expressed  in  the  form  of  geometric  coefficients 
which  have  the  same  general  characteristics  but  differ¬ 
ent  numerical  values  for  the  laminar  and  turbulent  cases. 
The  differences  result  from  different  values  of  two 
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characteristic  parameters  p  and  m  .  The  first  param¬ 
eter,  p  ,  controls  the  longitudinal  distribution  of  the 
ablation  intensity  and  appropriate  values  for  the  laminar 
case  are  -.5  and  for  the  turbulent  case,  -.2.  The  second 
parameter,  m  ,  controls  the  variation  of  the  ablation 
intensity  with  the  effective  local  deflection,  and  it  is 
shown  that  values  of  m  =  1  and  1.6  are  appropriate  for 
the  laminar  and  turbulent  cases,  respectively. 

An  arbitrary  but  constant  ^*'^se- shift  between  the 
oscillation  angle,  a  ,  and  the  corresponding  ablation 
effect  is  assumed.  The  determination  of  tne  appropriate 
phase  angle  in  a  particular  case  would  require  a  detailed 
analysis  involving  considerations  of  unsteady  boundary 
layer,  a  fluctuating  heat  input,  and  thermal  diffusivitv 
characteristics.  Such  an  analysis  was  not  attempted  and 
an  arbitrary  phase  lag  cf  20  degrees  was  assumed  for  the 
two  examples  used  to  illustrate  the  application  of  the 
method  to  the  laminar  and  turbulent  boundary  layer  condi¬ 
tions.  Some  characteristic  features  of  the  effects  of 
ablation  in  these  two  cases  are  pointed  out  and  the  re- 
suits  of  a  parametric  study  showing  the  individual  effects 
of  a  change  of  phase-lag  angle  or  oscillation  frequency 
are  ■ 1  lustrated. 
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SYMBOLS 

Area  of  cone  ''ase 

Proportionality  factor,  relating  the  flow 

deflection  produced  by  ablation,  3  ,  to 

the  mass  injection  parameter,  F  (Eq.  36) 

6 

P»u» 

Abbreviation  for  -pg— 

*e  e 

Coefficient  relating  the  local  rate  of 
convective  heat  transfer,  ^  ,  to  the 
local  angle,  A  and  x  coordinate 
(Eq.  29) 

Chapman-Rubes i n  parameter,  calculated  at 
the  reference  temperature 

Constant  coefficients,  determined  by  the 
cone  geometry  (Eqs.  11  and  12) 

Constant  coefficients  specified  by  the 
cone  and  ablation  geometry  (Eqs.  51, 

52,  53,  and  54) 

Constant  coefficients  determined  by  the 
cone  and  ablation  geometry  (Eqs.  91  and 
92) 

Local  skin  friction  coefficient 

Local  skin  friction  coefficient  without 
ab lat i on 

Fitching  moment  coefficient,  based  on  Ag 
and  d 

Static  pitching  moment  derivative,  per 
dC.JN 

rad  i  an(= 

Pitching  derivative  due  to 


q  ,  per  radian 
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Fe,  FW 


Lo'  LrL2 


ma,mc 


Pressure  coefficient  (= 


P  -  P. 


Specific  heat  (BTU/l b- °R ) 

Cone  base  diameter 

Nond imens i ona 1  mass  injection  parameters 
p  v  P  v 

given  by  p  y  •  and  p  y  ■  »  respective 


P  v 
rw  w 

Ke  e 


,  respectively 


Circular  frequency,  cps 
Effective  heat  of  ablation  (BTU/lb) 
Enthalpy  (BTU/lb) 

Abbreviation  for  integral 
Nond  imer.2  i  ona  1  mass  transfer  parameter 
(Eq.  75)  p 

Abbreviation  for  the  quantity  it 


ADDreviaiion  tot  cne  quantity  -n - ?rir“ 

Heffp»l» 

Nond imens i ona 1  factor,  appearing  in  Eq.  90, 
relating  the  local  skin  friction  coeffi¬ 
cient  to  the  mass  injection  parameter,  Fq 

Length  of  the  pointed  cone 

Axial  distances  from  the  cone  apex  to  the 
nose,  cone  shoulder  and  center  of  spheri¬ 
cal  cap,  respectively  (Figure  1) 

Mach  number  or  pitching  moment 

Rate  of  ablation  for  the  entire  body 
( lb/sec) 

Positive  real  quantity,  not  necessarily 
an  integer  used  as  an  exponent  to  the 
angle  A  (Eq.  29) 

Molecular  weight  of  air  and  coolant  gas, 
respect i ve ly 

Local  rate  of  ablation,  ifi  =  p,  v,  ( 1  d/ s q  ft-  sec) 

w  w  w 
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n 

P 


Q 

q 

* 

q 

Re 

r 

rN 

S 

s 

T 

t 

U 

v 

x,y,z 


X 


o 


x  i  *  *2 


P 


Integers,  used  in  binomial  series  expansion 
Static  pressure  or  exponent  expressing  the 
x-wise  variation  of  ablation  or  friction 
effects 

Heat  flux  for  entire  body  (BTU/sec) 

Dynamic  pressure  (=  Jr  ,  psf,  or 

angular  velocity,  radian  per  second 
Local  rate  of  convective  heat  transfer 
(BTU/sq  ft-sec) 

Reynolds'  number 

Radius  or  cylindrical  polar  coordinate 
Radius  of  spherical  nose 
Surface  area 

Distance  along  conical  surface 
Period  of  oscillation,  seconds 
Time,  seconds 
Stream  velocity,  ft/sec 
Normal  velocity,  ft/sec 

Cartesian  coordinates  with  respect  to  body 
axes  with  origin  at  the  cone  apex 
Distance  from  the  origin  to  the  axis  of 
osc illation 

Coordinates  of  the  leading  and  trailing 
edges  of  the  ablating  region,  respectively 

a  s  i  nqp 

Abbreviation  for  1 - t - 


Angle  of  attack,  radians 
Initial  angular  displacement,  radians 
Effective  angle  of  attack  to  account  for 
pitching  velocity  (=  a  +  t 

radians 

Local  flow  deflection  induced  by  ablation, 
rad ' ans 
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m 


* 


3  £ 

\  y 


Y 


6* 

K 


Tj 

6 

7k 

§ 

P 

E 

CT 

(P 

U) 


Ratio  of  specific  hects 
Net  local  deflection  from  the  free  stream 
direction,  for  the  oscillating  cone, 
rad  i  ar.s 

Boundary  layer  displacement  thickness 

Nond imens iona  1  factor  relating  the  change 

in  local  skin  friction  coefficient  to  the 

mass  transfer  parameter  (Eq.  78) 

Pressure  correction  factor  co  improve  the 

Newtonian  value  (Eq.  2) 

Cone  half-angle,  radians 

Phase-shift  angle,  radians 

Pluntness  parameter  (_  N'\ 

\  rB/ 

Fluid  mass  density,  slugs  per  cubic  foot 
For  s.  lation 
A 1 1 i tuao  dens i ty  rat i o 
Cylindrical  polar  coordinate 
Angular  frequency  of  oscillation,  radians 
per  second 


Subscripts 


a  Air 

9  Refers  to  cone  base 

c  Refers  to  coolant  gas 

e  «t  edge  of  boundary  layer 

l  Laminar 

N  Refers  to  the  spherical  nose  of  the  cone 

o  Without  ablation 

q  Indicates  partial  differentiation  wit!. 

respect  to  3t 
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s 

t 

w 

a 


os 


Stagnat i on 
Turbu 1 ent 

Conditions  at  the  wall  surface 
Indicates  partial  differentiation  with 
respect  to  .'ngle  a 
Free  stream  conditions 


Superset i pts 

(  )'  Value  for  the  pointed  cone 

*  VaUie  calculated  at  the  reference  temperature 


Additional  Symbols 


-(c) 

-(f) 

—  (F.P.  ) 

-U) 

—  (N) 
~(o) 

—  (t) 
~(P) 


Contribution  of  the  cone  surface 
Contribution  due  to  friction 
Flat  plate  case 
Laminar  case 

Contribution  of  the  cone's  spherical  nose 
Value  without  ablation 
Turbulent  case 

Contribution  due  to  ablation 
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APPENDIX  B 


DERIVATION  OF  SIMPLIFIED  EXPRESSIONS 
FOR  THE 

RATES  OF  HEAT  TRANSFER  TO  A  FLAT  PLATE 

AT 

SMALL  ANGLES  OF  ATTACK 


FluiDyne  Engineering  Corporation 


1 .  Laminar  Boundary  Layer 


The  expression 


s i na(cosa) 

T 


.5 


[.  186  +  .5  -jp  +  . 3 14 ( s i  not)2 


TT5 


(B 1 ) 


which  appears  in  Eq.  (19)  may  be  simplified  for  small 
values  of  a  by  using  the  expansio 

a3 

s  i  na  =  a  -  -s-T  +  •  •  • 

O  • 


and 


,  a2  a4 

cosa  «  1  -  -jj-j-  +  -fj  -  •  •  • 


and  applying  the  binomial  theorem,  neglecting  terms  in 
a4  as  fol  lows: 

/  J/2  /  a2  \1/2  12 

(cos (a) J  =  v1  -  JT  •••)  -  1  -  i  a  •" 


Then,  the  expression  B1  yields 


B 1 


(a  "  7l)  ( 1  "  T") 


[.  186  +  .5  ~  +  .314(a  -  yj)2] 


W 
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which  appears  in  Eq.  (24)  can  be  simplified  by  a  similar 
procedure  as  follows: 


[..97  +  .5^  +  -303 («-fr)  ] 

9 


a1,6(l  -  a2)  (l  -  .4  a2) 


T T 


(.  197  +  .5  -pp)  [l  + 


.303  a 


.  197  +  .5  -7T- 


#)] 


77 


3 

If,  as  before,  one  neglects  terms  of  order  fx 
then 


or  higher. 
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skin  friction  change  which  are  derived  from  flat  plate  data  modified 
to  account  for  the  conical  flow  conditions.  The  primary  independent 
quantities  other  than  the  cone  angle,  nose  bluntness  and  C.G.  locatio 
are  the  ablation  intensity,  geometry  and  phase-shift  w.r.t  the  oscil¬ 
lation  angle.  The  ablation  geometry  is  conveniently  expressed  by 
geometric  coefficients  of  the  same  general  type  but  different  numeri¬ 
cal  values  for  the  laminar  and  turbulent  cases. 

Two  examples  illustrate  the  application  of  the  analysis  to  the  lami 
nar  and  turbulent  boundary  layer  cases.  The  effect  of  ablation  was 
very  small  on  the  static  derivative  but  quite  considerable  on  the  dy¬ 
namic  derivative.  The  analysis  points  to  significant  differences  in 
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cases  with  potentially  serious  practical  consequences. 
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